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We analyze the dynamical chiral symmetry breaking in gauge theories solely by 
the non-perturbative renormalization group, without recourse to the Schwinger- 
Dyson method. First, we briefly review the basic notions and formulation, and 
clarify its great feature that it gives a systematic approximation scheme without 
any divergent series nor serious gauge dependence, compared to other perturbative 
and non-perturbative methods. Then we apply this new method to QED and QCD 
to find that our lowest level approximation improves the ladder Schwinger-Dyson 
equation results in a gauge independent way. We get the chiral phase structures 
and the critical exponents in the fixed gauge coupling analysis, and also calculate 
the chiral condensates in case of running QCD. 


1 Introduction 

Our aim is to attack the dynamical chiral symmetry breaking solely by the 
non-perturbative renormalization group (NPRG). The ladder Scliwinger-Dyson 
(SD) equation has been the best analytic toollo investigate it^a but is unsat¬ 
isfactory due to the strong gauge dependence^ and difficulty of going beyond 
the ladderp while the improved ladder approximation with the ruling gauge 
coupling constant gives good results in QCD even quantitativelycfQ However 
the improved ladder is just a model, having no firm theoretical base, thus there 
has been no way to further improve it. 

As everybody knows, NPRG must best fit the issue. We first define a 
renormalization group equation for some appropriate sub-theory space, then 
we get flows, critical surfaces and fixed points, and we can calculate critical 
exponents or anomalous dimensions etc. However nobody had known a good 
way of approximation to perform this strategy. Recently we found a good ap¬ 
proximation, which is a systematic and consistent method, contains the ladder 
and even the improved ladder SD results, thus we can go beyond the ladder, 
can improve the improved ladder, in a gauge independent way. It may be 
applicable to wide range of models and phenomena including supersymmetry, 
dynamical gauge symmetry breaking, topological effects, etc. 

“Invited talk at International Workshop on Perspectives of Strong Coupling Gauge The¬ 
ories (SCGT96), Nagoya, 1996. This talk is based on the work in collaboration with 
K. Morikawa, K. Shimizu, W. Souma, J.-I. Sumi, M. Taniguchi, H. Terao and M. Tomoyose. 
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In a battle field map of the renormalization group approaches (Fig.^, 
NPRG is just in the mid of the field, between the perturbative and the non- 
perturbative head quarters. It shares good features of both sides. The lattice 
simulation is in fact the most prominent in the non-perturbative regime, but it 
suffers serious problems of chiral fermions, and that it has exhausted the com¬ 
puter resources on the earth, although recently improved actions enhanced 
it even thousand times. It is not feasible to simulate systems with large hi¬ 
erarchies. Its great success in QCD is due to the fact that the confinement 
transition can be described within a relatively short range of energy scales. 
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On the other hand, in NPRG, there are no difficulty of chiral fermions nor 
large hierarchies. The most serious problem in NPRG is that it cannot respect 
the manifest gauge symmetry due to the introduction of the momentum cutoff. 
However this is not an essential difficulty, but just a practical issue. One should 
remember that the renormalization calculation in gauge theories could be done 
even with the momentum cutoff regularization, as was done for QED before 
the gauge invariant regularization was invented. 

Figure ^ shows a more personal view of the field. We, the Wilsonian 
people, is surrounded by other old countries. In other words we share all 
beauties of these countries simultaneously. We are now lonely attacking the 
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weakest neighbor of the republic of ladder also known as Nagoya, and finally 
will merge it, we hope. 

2 Wilsonian renormalization group 

In short, NPRG is a method to evaluate integral by solving a corresponding 
differential equation, 

F{t) = f{x)dx = n0) = 0, 

which is nothing but a definition of integral, though. We evaluate a path 
integral by snlving a functional differential equation which is called NPRG 
equationB’Dollj The essential ingredient is the Wilsonian effective action Seg 
defined by 

J Pcjiexp (-S'[(/)])= y exp (-S[(^>, (/)<]) = J !?(/)< exp (-Seff[(^<]), 

exp(-Seff[0<]) = y 2?0>exp(-S[0>,(^<]) . 

The high frequency modes are integrated out, and the effective action for 
the low frequency modes (^<) is obtained. This is equivalent to the so-called 
block spin transformation, and it physically describes the coarse graining or 
the transition from micro to macro. 

Now, we define the flow of SeS by continuously decreasing the cutoff A 
which is the boundary between the high and low frequency modes. The ef¬ 
fective action at A, Seff[A], is obtained after the high frequency modes (> A) 

have been integrated out: Z = 2?(/>exp (—S[(/); A]) , where the measure 

d 

contains low modes only. We evaluate the derivative of Seff[A], — Seff[0;A], 

where the fictitious time parameter t is introduced to control the cutoff: A(t) = 
exp(—f)A, i5A = A6t. Suppose we lower the cutoff by i5A, and evaluate 
the change of the effective action SSeS, which comes from the shell mode 
([A — (5A,A]) path integral. We expand SeS with respect to the shell mode 
field 0s, 

exp{-5'eff[0<; A - (5A]} = J T>(j)sexp{-Ses[(l)< -|-0s;A]} 

= exp(-5'eff[0<; A]) J l?0sexp|-y 
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where the momentum integration is carried over the shell momentum denoted 
by the prime. Only the first order term in must be kept to evaluate the 
derivative. The terms with more than two shell modes do not contribute to 
the first order, since they must be accompanied by at least two loop integral of 
the shell momentum, thus leaving 0((5A^). Then the shell mode path integral 
can be exactly carried out, which is just the Gaussian integral, resulting 




If [ f S^SeS ] SSeS ( \ ' 6S,S \ 


Graphically these contributions are represented by the following diagrams, ring 
and dumbbell, and its physical reasoning is clear. The shell mode path integral 
produces all the quantum corrections due to the shell modes, with external 
fields of the low modes. Note that the tree and the one loop diagrams complete 
the derivative, and no external momenta may flow in the loop except for the 
dumbbell endpoint vertices. 



Shell modes 
- Low modes 


To make NPRG equation physically sensible, we transform variables to be 
dimensionless, where the dimensional unit is taken to be the cutoff at that 
‘time’. This transformation runs the coupling constants, which is called the 
canonical scaling. Also we introduce the wave function renormalization to 
normalize the kinetic terms of every fields, which chooses a representative of 
the equivalent classes in the theory space. Then we finally get 
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where d is the space-time dimension, rj is the janomalous dimension. This 
equation, called the Wegner-Houghton equation^ describes the change of the 
effective action, that is, of a set of infinite number of coupling constants. Thus 
the right-hand side should be called a /3 functional, which controls the flows in 
the infinite dimensional theory space. Important is that it is obtained as an ex¬ 
act and closed form, therefore it is sometimes called the exact renormalization 
group. 


3 Approximation method 


Even though we obtained the exact renormalization group, it is impossible to 
solve such an infinite dimensional functional differential equation. Thus what 
is more important is if there is a good way of systematic approximation or not. 
The approximation method for the exact renormalization group is a projection 
to a sub-theory space. For example, we take a sub-theory space spanned only 
by the local potential terms in addition to the kinetic terms 


which is regarded as the lowest order of the expansion of the effective action 
with respect to the number of field derivatives. This is called the local potential 
approximation (LPA)pll3 where the corrections to the kinetic terms are ignored 
either, thus it gives vanishing anomalous dimensions. The Wegner-Houghton 
equation in LPA is reduced to read 


dt 


^[lnil + V'')]+d-V+^^V' , 


where V denotes the derivative with respect to the field (j), and Ad is a constant 
of the d-dimensional angular integral. We see the original functional differential 
equation is simplified to a 2-dimensional (t, (j)) partial differential equation, 
which is non-linear, and is to be solved as an initial value problem. Note 
that LPA still deals with an infinite dimensional theory space spanned by an 
arbitrary function V{4>). 

It is still not easy to analyze the system precisely. We reduce it to a 
finite dimensional sub-theory space by expanding the potential in terms of the 
polynomials in ip. For example, take a scalar theory with Z 2 symmetry in 3 
space-time dimension, and expand the potential as 


V{(l),t) = ai{t)p + a 2 {t)p'^ + a 3 {t)p^ ■ ■ ■ , p = 0^/2 . 
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Then the (3 functions (flow equations for the coeflicient functions) are given 
by 
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Even in this crude approximation with only 3 operators, we find a non-trivial 
fixed point realizing the ferromagnetic phase transition with the critical expo¬ 
nent V = 0.586, where the lattice simulation gives 0.629, such a good coinci¬ 
dence. We increase the number of operators or the dimension of the sub-theory 
space to see if the results converge. If it does, we understand we get the results 
for LPA sub-theory space. 

The approximation here is a projection to a sub-theory space, giving a 
projected /3 functions, which define a projected flow. The projected flow is 
different in general from the projection of the true flow. The true flow al¬ 
ways goes out of the sub-theory space. 

The projected flow pulls it back on the 
sub-theory space in each step of trans¬ 
formation. In some exceptional case as 
the large-iV theory, the projected flow 
coincides with the projection of the true 
flow if we choose a good coordinate sys¬ 
tem of the theory spaca, which we call 
the perfect coordinatell3 

Generally we may expect that enlarging the sub-theory space, the physical 
results converge, or oscillate at worst. This is a natural expectation. Enlarging 
the sub-theory space is not directly related to any divergent series expansion. 
It is more like enlarging the total lattice size in simulation, where nobody 
expects the results to diverge. However, usually we must think out a good 
series of sub-theory spaces to get fast convergence. We have learned that an 
‘environmentally friendly’ coordinate system assures better convergence. For 
example, expanding the potential, we might de:^e the following two schemes: 
a fixed scheme (A) and a comoving scheme (B)PJ 


Projected Flow and True Component 



A : V{p) = Y, 

n—1 


B: H(p)= ^ a„{p-po(i))” 

n—1 
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The comoving scheme is a dynamic coordinate since its origin of expansion po{t) 
is taken as a minimum of the potential at time t. Physically it is better since it 
uses more appropriate prop^a tors a t and it actually improves convergence 
drastically as is seen in Fig.^J’U’E3[| 

4 Comparison with other non-perturbative methods 

We compare our NPRG results with other standard non-perturbative methods. 
First let us examine the leading cxijtical exponent of d-dimensional scalar theory 
comparing with the e-expansionll3 As is seen in Fig.^, the e-expansion gives a 
divergent series. Of course it is extensively investigated that these divergent 
series can be Borel resummed to give a convergent value which is also displayed 
in the figure. Our NPRG analysis gives a globally right behavior in the lowest 
order approximation LPA. 


Truncation dependence of the leading exponent 
D=3 single scalar, Evaluated by LPA W-H eqn. 
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Next we analyze scalar theories with N-components in 3 .space-time di¬ 
mension, and compare NPRG with the 1/A^ expansion (Fig.^)E3ll3’0 The 1/A^ 
leading result is just v = I, and we show the next to leading (0(1/A^)) and 
the next next to leading {0{1/N'^)) results. It is clearly seen that the 1/A^ 
expansion is also a divergent series. For small N theories, the leading result 
(i/ = 1 classical result) is the best, and adding higher orders in 1/A^ makes it 


After the conference, we find that even the comoving frame gives finally diverging results. 
However, it is at extremely large or^rs, thus it is after the highly accurate convergence of 
even 8 digits for u has been obtaineco. We thank T. R. Morris to motivate us to investigate 
very large orders. 
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worse. Our NPRG analysis gives a stable and globally correct results also in 
this comparison, and becomes exact at —> oo limit. We expect that when 
enlarging the sub-theory space beyond LPA, we get closer to the true values. 
Some of the beyond LPA results are also plotted in Fig.^ and our expectation 
proves right at least in the small N region. Note that Borel resummation for 
1/A^ expansion is not feasible at all. 


Comparison with 1/N expansion 



What is a crucial difference between NPRG and, for example, 1/N expan¬ 
sion? The 1/A^ expansion contains only limited structures of diagrams corre¬ 
sponding to the order of the expansion, while NPRG contains every structured 
diagrams, exactly for the leading, but not exactly for higher orders in 1/N, 
even at the lowest level of LPA. Generally NPRG contains every powers of 
parameters as e, 1/A^, and coupling constants, and it contains the lowest order 
exactly. This feature assures that NPRG gives globally correct behavior, and 
it is expected to converge when enlarging the sub-theory space. 

We show another striking results in case of sca]^ theories in 1 space-time 
dimension also known as the unharmonic oscillatorPJ whose action is given by 


S = 


dt 


± -x^ + Xx 
2 2 



We evaluate the energy gap between the ground state and the first excited 
state using the infrared effective potential obtained by NPRG. In case of the 
single-well potential (Fig.^, our results are almost exact for whole range of the 










coupling constant A, whereas the standard perturbation series shows its diver¬ 
gent nature even at the weak coupling region (although Borel resummation is 
possible in this case). 


Effective mass of single-well anharmonic oscillator 
W-H RG Equation : cutoff scale = 100 




Fig.| 


Fig.| 


In case of the double-well potential, we compare our results with the instanton 
calculation (Fig.||). At A ~ 0.08 the first excited state energy crosses down the 
central barrier, where the LPA results start deviating from the exact values. 
At this point the instanton results are already unreliable. The LPA works 
extremely well except for the deep region where the tunneling rate is very 
small and the dilute instantons become exact. 

To see how NPRG works well for the double-well potential is very im¬ 
portant, since it is related to a more general issue of how NPRG can treat 
essential singularities and topological objects in field theory. The LPA gives 
a convex effective potential in case of the spontaneous symmetry breaking, 
which indicates that LPA is able to take account of the_global domain wall 
structures without which the convexity may not appearlld However it seems 
that LPA is not successful to evaluate precisely the action accompanying such 
topological configurations since it fails to give the correct energy gap in case 
of low tunneling rate. 


5 Non-perturbative RGE vs. Perturbative RGE 

Here we recapitulate the features of NPRG, compared to the standard per¬ 
turbative RGE (PRGE). PRGE defines the /3 functions which are calculated 
by the loop expansion, giving the leading log series expansion or the improved 
perturbation, while NPRG (3 functional is infinite dimensional, but it is exactly 
given by tree and one-loop diagrams. Multi-loop effects come out through the 
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higher dimensional operators. PRGE gives the asymptotic series at best, and 
actually other non-perturbative methods including the e-expansion and the 
1/A^-expansion do as well, while in NPRG the approximation is defined by 
a projection to a sub-theory space and systematic enlargement of it, thus no 
diverging series may appear since no expansion in any parameter is employed. 

The standard effective action r{(j)) obtained by PRGE does not serve the 
dynamical symmetry breaking due to a composite order parameter. We need 
to evaluate Cornwall-Jackiw-Tomboulis (CJT) effective potential for the conn 
posite operators or equivalently to solve the Schwinger-Dyson (SD) equation!] 
However, NPRG automatically contains the CJT effective potential through 
the Wilsonian effective action including higher dimensional multi composite 
operators, thus the non-perturbative vacuum can be directly analyzed. 

NPRG cannot manifest the gauge invariance due to the momentum cutoff. 
Then theory space must span gauge non-invariant operators as well, and we 
have to pick up the gauge-upariant theory space (which differs from the gauge 
invariant operator space) c2Hl3 This is a serious problem for load of the practical 
calculation, but it is not an essential difficulty and can be treated anyway. 

To summarize, NPRG method is very general, robust and global, and no 
divergent series appears. Recent development has been done by numerical 
integration of the NPRG equation in LPA and beyond LPA sub-theory spaces. 
This is a balanced research between analytic and numeric. The lowest order 
sub-theory space, LPA, has exhibited the features that it contains one-loop 
perturbative results, exact in the large N leading, exact in the 1st order of 
e-expansion, and furthermore contains the improved ladder SD results. 

6 Dynamical chiral symmetry breaking 

Due to shp£.tage of pages I have, we just briefly describe the basic results of 
our worksEj We consider a chiral invariant gauge theory and set up NPRG 
equations in LPA. We include multi-fermi operators, 4-fermi, 8-fermi, • • •, re¬ 
specting the chiral invariance. The (3 functions for 4-fermi operators do not 
depend on 8- or more-fermi operators due to the chiral invariance. Thus our 
multi-fermi coordinates is a perfect coordinates and the flow of 4-fermi oper¬ 
ators can be solved exactly in LPA, coupled with the running of the gauge 
coupling constant e (A = 3e^/47r^): 
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Figure ^ shows the diagrams contributing to the (3 functions of the scalar 4- 
fermi Gg and vector 4-fermi Gy interactions, where the dashed lines represent 
the gauge boson. 



Fig.| 


P-Function of 4-Ferini Coupling 



Fig-0 


We define a part of the (3 function, the ladder part, as in Fig||. If one takes 
only this part in the following analyses, one obtains exactly the same results 
as those by the ladder SD equation in the Landau gauge. (This uniqueness of 
the Landau gauge comes from the fact that the fermion anomalous dimension 
vanishes in this gauge.) Remember that the ladder SD suffers a strong gauge 
dependence^ whose reason is clear here. The ladder picks up just a part of the 
gauge invariant set of diagrams in the (3 function, taking the ladder but not the 
crossed ladder. On the other hand, NPRG /3 function takes both, no difference 
for either ladder or crossed ladder, thus it respects the gauge independence of 
flows and of the physical results. 

One readily sees (Fig.|l^) that at A = 0, the scalar 4-fermi flows has an 
ultraviolet fixed point which is nothing but the Nambu-Jona-Lasinio (NJL) 
criticality, and an infrared fixed point of the Gaussian one. There appear two 
phases, strong and weak. Note that the 8-fermi operator has no fixed point 
without gauge interactions, indicating the triviality of the pure fermi theory. 
When the gauge coupling is switched on, the Gaussian fixed point moves up, 
the NJL fixed point moves down, and they finally meet together to pair an¬ 
nihilate at a critical gauge coupling Ac, above which there is only one phase 
left. To identify these two phases, we investigate the CJT effective potential 
for the chiral condensate < iliip >, and find that the weak phase is a symmetric 
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phase, while the strong phase exhibits the spontaneous symmetry breaking of 
the chiral symmetry, where the potential takes the flat bottomed convex shape 
characteristic to the spontaneous symmetry breakdown. This structure is dis¬ 
played as a phase diagram in the scalar 4-fermi vs. the gauge coupling plans 
in Fig.^ which should be compared with the old traditional phase diagrama 
We may calculate a physical quantity of this system, the anomalous dimension 
of the mass operator '0'!/) (Figjl^). It is enhanced compared to the ladder SD 
results. Note that our results are LPA exact, including ‘beyond the ladder’ 
diagrams, and are gauge independent. 


Flow of 4-Fermi Coupling 


Anomalous Dimension of 


Case of Ladder part 



Fig0 


Fixed Gauge Coupling 



Fig.0 


We may run the gauge coupling constant. The U(l) gauge interactions 
keep the two phase structure. However the fixed points are NJL and the 
Gaussian only, and the renormalized trajectory_is on the trivial line, which 
is consistent with the recent lattice simulationsE^ Asymptotically free gauge 
interactions wipe out the symmetric phase, and only the Gaussian fixed point 
survives. Therefore the whole region of the theory space belongs to the chiral 
symmetry broken phase. We can evaluate physical quantities, the chiral con¬ 
densate and the fermion mass. In this running gauge coupling case, we proved 
that the ladder part /3 function now gives exactly the same results as those 
by the improved ladder SD equation in the Landau gauge, a la HigashijimaP 
where the scale parameter of the running gauge coupling constant takes a larger 
momentum of the fermion legs. This particular definition of the running gauge 
coupling constant exactly emerges from the NPRG LPA with the ladder part. 
Now we know what is included and what is missing in the improved ladder SD 
equations, thus establishing the improved ladder for the first time. 

To evaluate the infrared physical quantities, the bare multi-fermi coor- 
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dinate system is not good, since it uses the massless fermion propagators. 
We introduce a collective coordinate of a scalar field coupled to the 'ip'4’ op¬ 
erator, which is another example of the environmentally friendly coordinate 
system. Within the ladder part, we have calculated the chiral condensate and 
the fermion mass as in Figjl^ which shows good convergence within low dk 
mensional sub-theory spaces, and agrees with the improved ladder SD resultsQ 
Taking the non-ladder /3 components into account, we proceed to evaluate in¬ 
frared physical quantities beyond the ladder SD results, in a gauge independent 
way. That will be the first results of gauge independent and including ‘beyond 
the ladder’, thought it is not an easy task. 


Dynamical Chiral Symmetry Breaking in QCD 

Dynamical mass 1(0) and Order parameter<^> 



Fig.y 


To deal with the gauge interacticpsL_more seriously, we need to formulate 
NPRG with a smooth cutoff functionEJ’M Then the cutoff scheme (the profile 
of the cutoff function) dependence might cause a problem. However we find 
that indeed such scheme dependence affects phase criticalities (the critical cou¬ 
pling constants) which are themselves not physical quantities, but the physical 
quantities like the anom^us dimension of the mass operator has only very 
small scheme dependenceo 

We have clarified that NPRG equipped with the sub-theory space approx¬ 
imation potentially solves the dynamical chiral symmetry breaking in gauge 
theories, which can go beyond the ladder SD results restoring the gauge in¬ 
dependence. This new method of treating non-perturbative QGD will give us 
a promising tool alternative to the lattice simulations. However we may need 
some revolutionary new techniques of evaluating the /3 functional to solve QCD 
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with larger dimensional sub-theory spaces of the gauge sector keeping the gauge 
invariance. 
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